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Abstract 

We apply the recently proposed transfer matrix formalism to 2-dimensional quantum 
gravity coupled to (2, 2k — 1) minimal models. We find that the propagation of a parent 
universe in geodesic (Euclidean) time is accompanied by continual emission of baby universes 
and derive a distribution function describing their sizes. The k —>■ oo {c —>■ — oo) limit is 
generally thought to correspond to classical geometry, and we indeed find a classical peak 
in the universe distribution function. However, we also observe dramatic quantum effects 
associated with baby universes at finite length scales. 
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1 Introduction 



One fascinating aspect of quantum gravity is the possibility of topology changing processes 
where a compact connected 3-geometry (a "baby universe") splits off or joins a large "parent 
universe." Such processes, peculiar as they may seem, are quite natural at the Planckian 
distance scales, where geometry undergoes large quantum fluctuations. Even if the topology 
changing processes are confined to the smallest distance scales, they may drastically affect the 
observed physics. As shown in ref. , summation over all possible emissions and absorptions 
of the Planck-size baby universes effectively gives rise to averaging over the fundamental 
constants of nature. If the weight in the average is sharply peaked, then the fundamental 
constants are determined 0. The wormhole theory of fundamental constants shows that, via 
some unexpected connections, the short-distance properties of quantum gravity can affect 
the universe at large in crucial ways. The status of this theory when applied to 4-dimensional 
quantum gravity is somewhat uncertain due to the instability of the Euclidean path integral. 
It may be, for instance, that the weight in the average over the coupling constants is not 
sharply peaked 0. If, on the other hand, the weight is sharply peaked, then it is not 
quite clear what prevents the emission and absorption of macroscopic baby universes [Q. 
While these puzzles await their resolution through a better understanding of 4-dimensional 
quantum gravity, we may try to gain some intuition about the topology changing processes 
from low- dimensional models. 

The 2-dimensional Euclidean quantum gravity coupled to c < 1 conformal matter appears 
to provide some promising toy models which are not plagued by any perturbative instabilities. 
Our understanding of these theories has dramatically improved recently, in part due to the 
success of the matrix model techniques [Q. With their help we can perform summations over 
discretized random surfaces and extract the exact results in the universal continuum limit. 
Until recently, however, the available exact solutions have not shed any light on the physics 
of topology changing processes. The problem is that the matrix models do not give any 
direct information on the internal geometry of 2-dimensional space-times. Some first pieces 
of information were extracted via direct Monte Carlo studies of triangulated random surfaces 
0. It was noted that the set of points at a given geodesic distance D from some point P 
typically consists of many disconnected loops, whose number grows rapidly with D. The 
study in ref. gave a first indication that, as a 1-dimensional universe (a string) propagates 
forward in geodesic (Euclidean) time, it is very likely to continually emit baby universes. This 
phenomenon received a detailed quantitative confirmation in a recent very interesting paper 
by Kawai, Kawamoto, Mogami and Watabiki |^. Relying on new combinatorial techniques 
that they invented, these authors derived a formula for the average number of loops with 
lengths between L and L + dL located at a geodesic distance D from some point on the 
surface, 

p(L, D)dL = — 1^ (x-^/^ + la:-^/^ + ^x^^) e^rfL (1) 
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where x = L/D"^. A number of beautiful conclusions follow from this formula. First of all, 
p{L, D) is not normalizable at small L, which means that the microscopic baby universes are 
overwhelmingly likely to be emitted. Secondly, apart from an overall dimensionful factor, p 
is a function of the dimensionless scaling variable x. Introducing the moments, 

{^0= / p{L,D)L"dL, (2) 
Jo 

one finds that for n < 2 they are dominated by the non-universal short- distance cut-off. But 
for n > 2 there are simple scaling relations {L^)d ~ D'^^, whose form follows essentially 
from D having the dimension of (Length) 

The distribution function p{L, D) is very important because it quantifies the effects of 
baby universes in pure Euclidean 2-dimensional quantum gravity. The next question is how 
the coupling to matter affects the topology changing processes. In this paper we begin 
to address this question by considering the emission of baby universes in 2-dimensional 
quantum gravity coupled to the {2,2k — 1) minimal models, whose central charges are c = 
1 — 3(2A; — 3)^/ (2A; — 1). These theories have been identified [§] with the multicritical points of 
the one- matrix model 0,0. The k = 2 theory corresponds to pure gravity, where it suffices 
to consider the discretizations of random surfaces with squares only. For = 3 (c = —22/5) 
one needs squares and hexagons; for /c = 4 one needs squares, hexagons and octagons; etc. 
For k > 3 some polygons enter with negative weights, which is not surprising given that 
these models are not unitary. We find that, for odd k, the negative weights have so much 
effect that there is no sensible positive p{L,D). For even k, however, there does exist a 
positive p{L,D), whose calculation constitutes the main result of this paper. We find that 
it depends on the scaling variable x = L/ D^^^^^^/'^^ and for small x diverges as x^^^^^^. 
Thus, the emission of microscopic baby universes becomes more enhanced with increasing 
k. For large fc, c ^ —oo and the sum over surfaces is expected to be dominated by classical 
geometry. In this limit we indeed find that p(L, D) exhibits a sharp peak that corresponds 
to the macroscopic classical geometry. Surprisingly, we also find that baby universes remain 
prevalent for all length scales less than a fixed constant, Sc, times the macroscopic (classical) 
length scale, and calculate the critical value Sc- As a further application of these ideas, we 
discuss the production of baby universes by a very large parent universe of length Lq. We 
find that their distribution is governed by a very simple scaling law, 

p(Lo ^ oo, L, D) ~ LoDL~'-'/^ , (3) 

which applies to both even and odd k. 

The structure of our paper is as follows. In section ^ we review and rederive the neces- 
sary matrix model results. In section ^ we extend the transfer matrix formalism of ref. 
to general discretizations, which are necessary to describe the theories with matter. In sec- 
tion ^ we calculate the transfer matrix in some limiting cases, and in section ^ we study the 
branching structure of space-time by deriving the distribution functions for baby universes. 
In section ^ we conclude with a few remarks. 
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2 The Disk Amplitudes 



For the calculation of the transfer matrix, the only result needed from the matrix models is 
the disk amplitude — that is, the partition function for surfaces with one boundary loop. It 
was shown in that in the k*'^ multicritical theory, the universal part of the disk amplitude 
is 

r) = i(v^)^-V2^,_i/2(v^L) (4) 

where r is the cosmological constant and Kk_i/2 is a modified Bessel function. We will 
need to keep track of certain non-singular parts of the disk amplitudes because they will be 
important in the calculation of the transfer matrix. Thus, we present our own calculation of 
the disk amplitudes. 

Throughout this paper, we will be working with surfaces without handles, so it is sufficient 
to use saddle point techniques in the matrix models. The partition function of the k^^ 
multicritical model (that is, the sum of the weights of all vacuum diagrams, including the 
disconnected ones) is 



Z = J exp [-(3 Tr ^ ($)] (5) 
where $ is a x matrix with ~ /5 — > oo, and 

vi^) = E ^^'"^ (6) 

m>l -^''^ 

is a potential energy function. (We restrict ourselves to even V{^) because they are more 
easily analyzed, but we believe that more general V{^) do not exhibit any more general 
behavior.) 

Let u{X) be the density of eigenvalues of normalized so that / dXu{X) = 1. As shown 

in 



121 , u{X) has support [—y/z, ^i] for some z, and 



For A G C — [— i/i, i/i], define 



N fV^ , u(u) 



F(A) is analytic except for a branch cut along [— -^/i, -y/i]: for A G [— ^i, y^] and e > 
infinitesimally small, F{X ± te) = V'{X) T 2m{N/p)u{X). Also, F(A) = 2{N/l3)/X + 0{X-^) 
for large A. It is not hard to prove that these properties uniquely determine F{X) and hence 
u{X). It turns out that F{X) has the form 



F(A) = r'(A) - /(A)Ayi^ (9) 
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where /(A) is a polynomial. Both z and /(A) are determined by the large A expansion of 
F(X). Defining for convenience the functions g[X) = XV'{X) and h{X) = A^/(A), we find that 



>W (l - - (E c.„A-j (e §11;) ^ MA) + 2| . O(A-), (10) 

The 0(A°) terms determine z: 

y c ^z- = 2^ (11) 

We can express the partition function Z as a power series in the couplings Cm and 
the parameter N/(3. If we consider the potential V and hence the Cm to be fixed, then 
Z has a radius of convergence in N/(3\ we construct the potential so that this radius of 
convergence is 1. For N/(3 > 1 — that is, for (renormalized) temperatures higher than a 
critical temperature — the series expansion of Z diverges. The behavior of Z near the critical 
point is what describes the quantum geometry of random surfaces, since at the critical point 
the behavior of Z is dominated by Feynman graphs of large order. It turns out that the 
singular behavior of Z at the critical point is determined by the dependence of z on N/ (3 in 
Eq. ([TTD . The /c*^ multicritical model is constructed by finding Cm such that Eq. (0) takes 
the form 



for some Zc. We shall use 



On 



\m+i 2A;! m! 
^"^^ (2m)!(A;-m)!A;™ " ^ (13) 

for m > k 



which gives Zc = 4:k. 

Let e be the lattice spacing of our random lattice. Then area is measured in units of e^, 
which will be taken to zero as the average number of plaquettes in the discretizations of the 
random surface diverges, so as to yield finite area. The parameter /iq corresponds to the 
lowest dimension operator in the k^^ multicritical theory, and it was at first believed that 
if one wrote fiQ = {2e)^tQ, then to would be the cosmo logical constant. This is false for the 
k > 2 models, which are non-unitary: to has dimensions of (Length)"'^, and is the coupling 
constant for the gravitationally dressed conformal field of the lowest dimension. The true 
cosmological constant is by definition conjugate to area, so it has dimensions of (Length)^^. 

We must consider a more general perturbation around criticality. Specifically, we perturb 
the potential by replacing Cm — > Cm + (^2m)\km f^m for all m > 0. Eq. ([llD then reduces to 

i-f)'^E.,.(f)". m 
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Higher dimension scaling operators are introduced by choosing the /i^ so the that right 
hand side becomes tie''~\l — z/zcY for some / < k, where ti is the couphng to the l^^ 
scahng operator. As was shown in []T0|, the couphngs to the operators corresponding to 
conformal fields in the Liouville theory are analytical combinations of ti with a definite 
overall dimension. The cosmological constant r corresponds to a perturbation with ti = 
Q;/2'^~'r'''^~'''/^ for Z — even and < I < k — 2. The ai are real numbers chosen so that the 
singular part of the one-loop amplitude obeys a Bessel equation, which emerges from the 
Wheeler de Witt equation. 

(l-f)' = E«K2ev^r'(l-f)' , (15) 

and by requiring Z); ct/ = 1 (which is equivalent to fixing a normalization for r) we find 
z = Zc{l — For k = 2 and k = 3 there is only one ai, so r cx tk-2- For = 4, it is 

found that ao = —1/5 and a2 = 6/5. We do not know a general way to determine the ai 
except by straight calculation. 

The partition function for surfaces with one boundary is the continuum limit of the 
Green's function of the field theory above. The /-point Green's function is 

Gi = {Tt^^) = f^' dXu{X)XK (16) 
It is convenient to introduce G{y) = J2i>oGiy'', since it is easy to see that 

G{y) = J^Yy^ (l/y) = {9 (1/y) - h {l/y) - zy^) . (17) 

Since aX z = Zc the radius of convergence of G{y) is yc = Xj^pZc-, the continuum limit is 
taken by setting y = ycexp{—e(), where now ( is conjugate to the length of the boundary. 
The continuum limit (e — 0) of G{y) was calculated using Mathematica, with the following 
results: 



k = 2: GiO = I - |C6 + ^6^/2 /2(C, r) + 0(6^) 
where ^(C, r) = 2 (C - ^v^) + 



k = 3: GiO = f - tC6 + <^^e^ + '-^e'/'UC r) + 0{e^ 



where /3(C, r) = | (-C + Kv^ + i^) \/C + 



k = A: G(C) = f - + ^6(-y3C^) ,2 + 32C(17|-185C!) ,3 + 51^,7/2 f^^^^ ^) + Q' 



35^ ' 35 ' 525 '7 

l(e-\ev-r-\cr+\ 



e 



where ^(C, r) = f - K'v^ " + ^'^^ VC + 



(18) 

The coefficient fk{C, t) of the leading nonanalytic term in e is the universal part of the disk 
amplitude. /^(C, t) is, up to a numerical factor, the Laplace transform of the disk amplitude 
Fk{L, t) in Eq. (^). The lower order analytic terms in ( and r correspond to zero length and 
zero area terms, and most of them can be dropped, for the following reason. We are free 
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to adjust ( by an analytic function of ( and r, which, in order to avoid trivial additive and 
multiplicative rescalings of (, and to preserve the dimension of (, should have the form 



C ^ C + ner + r2eC + rgeVC + r^e^C + ■■■ (19) 

where the are c-numbers. Such a redefinition preserves the universal term fk{C,T) and 
corresponds simply to a different way of treating zero area and zero length terms, as pointed 
out in |]TU| . By an appropriate choice of the in Eq. (|TP|) , we can absorb all but the first 



two analytic terms of the expressions in Eq. (p!8| ) into (. We then find the simple form 

G(C) = ^ - 7) + "^'"'^'/^(C' ^) + ^(^') (20) 

where 



a = k-3/2 (21) 

and a is a numerical factor chosen for each k so that the leading term of a small r expansion 
of/,(C,r) is (-1)V+V^- 

For later convenience, we mention one more mathematical point: since z/zc = 1 + 0(e), 
we can alter the perturbation Eq. (^) to read 




where n is a fixed integer, and the equation z = Zc{l — 2eA/r) will receive corrections that 
are analytic in e and of order and higher. The non-universal terms in Eq. (|18D will change, 
but the leading analytic and leading nonanalytic terms will be unaffected, so by making an 
appropriate redefinition of the form Eq. (|19D, we still arrive at Eq. (^). If we take n = 2, 
then fiQ = and /ii = 0. This is desirable because it means that the weights of planar 
Feynman diagrams do not depend on the number of edges. 



3 The transfer matrix 

The focus of is the analysis of the evolution of a loop through some fixed geodesic distance 
on the surface, where on a discretized surface geodesic distance is defined as the minimal 
number of plaquettes one must traverse to get from one point to another. The goal is to 
calculate the partition function of a tube with one entrance loop and one exit loop, such 
that each point on the exit loop is a fixed geodesic distance D from the entrance loop. 

On a discretized surface, one starts with an entrance loop 7 (see Fig. 1), defines a "for- 
ward" direction for geodesic distance (inward in Fig. 1), and thinks of advancing the loop 
along the lattice one step at a time. To accomplish this "one-step deformation," as it was 
called in 0], one first removes any double links that may exist on 7 and then moves each 
remaining link across the plaquette it borders in the forward direction. Clearly, this process 
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will sometimes split 7 into several loops. To phrase it another way, once we have removed 
all double links from 7, we color all the plaquettes that the remaining links border in the 
forward direction. Each side of a colored plaquette that is not part of the entrance loop 
7 becomes a link in one of the exit loops. Since we are interested in having just one exit 
loop, we designate one of the loops evolved from 7 as the exit loop, and close off the other 
loops with disks generated by the disk amplitude of the matrix model. A series of one-step 
deformations gives an exit loop which is a fixed geodesic distance from the entrance loop. 

We now want to ask the following combinatorical question: given that the entrance loop 
7 has I links and the exit loop 7' has V links, how many ways are there to construct a 
discretization between 7 and 7', using a given number of squares, hexagons, octagons, etc., 
such that 7' is the exit loop evolved from 7 in the course of a one-step deformation, as defined 
above? Or, in terms of the Feynman diagram which is the dual graph of the discretization, 
how many planar diagrams can be drawn on the surface of a cylinder with / external legs 
pointing "down" and /' external legs pointing "up", with specified numbers of vertices of 
order four, six, eight, etc., and such that every "up" leg is connected to the same vertex 
as some "down" leg (but not necessarily vice versa)? To answer this question, we assign a 
multiplicative weight g^n to each 2m-gon (equivalently, to each vertex of order 2m on the 
dual lattice) , assign as an overall weight to each of the discretizations between 7 and 7' the 
product of the weights of the plaquettes used to build it, and let Ni^i^ be the sum of the 
weights of all such discretizations. Ni^ii is the generating function that answers the above 
combinatorical question, but of course its more interesting property is that it is the partition 
function of tubes, or rather ribbons, of geodesic width one (in lattice units) and with entrance 
and exit loops of / and /' links, respectively. 

We propose to evaluate 

N{y,y')^ E N,,ryyf (23) 
i,i'>i 

using combinatorics and the disk amplitude from the matrix models. For calculational 
convenience we will designate the entrance loop as unmarked and the exit loop as marked. 
On the dual lattice, this corresponds to considering the / external legs of the planar diagram 
that point "down" to be distinguishable only up to cyclic permutations and the /' external 
legs that point "up" to be completely distinguishable. A different convention of marking 
would change Ni^i> only by a factor of I or V . The Green's function G{y) of the matrix field 
theory refers to diagrams in which the external legs are distinguishable — that is, G{y) is the 
amplitude of a disk with a marked boundary. The weight Qm assigned to a 2m-gon by the 
Feynman rules for the matrix field theory is 

- ^^'^ 

Let us first consider the case where only squares are used in the discretizations. All the 
discretizations contributing to N{y, y') may be built up as follows (here we are paralleling 
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closely the work of only using squares instead of triangles). The marked point on the 



exit loop must border one of the shapes in list a) of Fig. 2: the curved loops on the last three 
shapes denote the insertion of a disk. As we follow the exit loop around its length, we can 
encounter any of the shapes in list b) in any order: the entrance and exit loops are connected 
from the corner of one shape to the next, and in the end the last shape is connected back to 
the original member of list a), closing the loop. 

It is perhaps not transparent why the weighting of each shape is what it is, so let us work 
through the example shown in Fig. 3. The loop above the square indicates the insertion of 
a disk amplitude, which must be marked because something must pick out the point where 
the disk boundary meets the square. The disk discretizations with less than two sides on the 
boundary are omitted because there have to be two boundary sides that match with the two 
free sides of the square. Let us think now in terms of the dual lattice. Given any Feynman 
graph contributing to the disk amplitude, we form a graph of the desired type by attaching to 
an adjacent pair of its external legs another vertex. The weight of the original disk amplitude 
graph must be multiphed by g2 to get the weight of the new graph, because we have added 
one vertex and tied up two external legs while adding two new external legs. Summing over 
all allowable disk amplitude graphs (those with at least two external legs) then gives a total 
weight of g2{G{y) — Gq — Giy), as claimed. Incidentally, the first, "undrawable" term g2y'^ 
in Fig. 3 corresponds to the two free sides of the square being identified. It might be helpful 
for the reader to identify which disk amplitude diagrams correspond to each term in Fig. 3. 

Now it is not hard to write down the weight that would be assigned to a 2m-gon bordered 
on n sides by the exit loop. Starting with a disk amplitude graph with at least 2m — n — 2 
external legs, we add to it one 2m- vertex, tying up 2m — n — 2 of its external legs and adding 
n + 2 new external legs, 2 of which are entrance legs and n of which are exit legs. Thus the 
total weight is 



Note that because of the planarity of the surface, the sides of the 2m-gon bordered by the 
exit loop must be contiguous. When 2m-gons are allowed in the discretization, we include 
all the 2m-gon shapes in list b), and also in list a) with a multiplicity determined by the 
number of exit links on the shape. 

The combinatorical problem is now quite simple: we construct an arbitrary discretization 
from one member of list a) and some sequence of members of list b). The total amplitude 
when only squares are used is 




(25) 



N{y,y') = {3y''yg2 + 2y'Y92G{y)+y'yg2{G{y)-Go)) 



X 



oo 



T.[y"y92 + y'V92Giy) + y'yg2{G{y) - Go) + g2{G{y) - Go - G,y) + y'G{y)] 



n 



n=0 
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^ ^y''yg2 + 2yWg2G{y) + y'yg2{G{y) - Gp) 

1 - y'^yg2 - y'Wg2G{y) - y'yg2{G{y) - Go) - g2{G{y) - Go - G,y) - y^G{y) ' ^ ^ 

Expanding out the first expression gives terms corresponding to individual discretizations. 
Let us define 5* to be 1 minus the sum of weights in list b): 



S = l- y'hg2 - y'Vg2G{y) - y'yg2{G{y) - Go) - g2{G{y) -Go- G.y) - y'G{y). (27) 

Recalling that each shape from list b) enters into list a) with a multiplicity equal to the 
number of exit links on the shape, we find that 

a logy 

In the general case where arbitrary 2m-gons are allowed, list a) and list b) are appropriately 
expanded, and 



C2m—1 / 2m— n— 3 \ \ 

y''^-'ygm+j:y"'y''''"-""'9m{G{y)- E Gw']] 

2m-2 / 2m-n-3 \ 

= 1 - yy' + -,g{y') - 2y'G{y) -y'T^^Y. (y'yT G{y) - G,y' ' (29) 

y m>l y n=0 \ 1=0 / 

and Eq. (|^) still holds. 

The continuum limit of N{y, y') is taken by expanding about the radius of convergence 
of each of the variables. The convergence of 5* is determined by the convergence of G{y), 
which we explained in Section ||: ?/c = and we set y = ycexp{—e(). For fixed y < yc, 

S is entire in y', but N{y, y') is analytic in y' only up to the magnitude of the zero of S 
nearest = 0. It turns out that for C = and r = 0, this zero is at y' = . Hence we set 
y' = exp(— eC). An analytic redefinition of of the form 

/2 



C-^C + rier + r^eC' + r^eCC + r^^C + • • • (30) 



is allowed, for the same reasons as for the redefinition Eq. (^9|). Thus when we expand 
N{y, y') in e, we need only retain the leading analytic and leading nonanalytic terms in 
C'. Writing N{CX') place of iV (?/cexp(— eC), ?/(r^ exp(— eC')); we find remarkably simple 
results: 

where a' is another numerical factor. 
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Let Ni^i'{d) be the lattice partition function of tubes of geodesic length d with an un- 
marked entrance loop of / links and a marked exit loop of /' links. The great insight of is 
that this object has a simple composition law: 



oo 



NiAdi + d2) = E Ni,i,idi)Ni.,,{d2) . 



(32) 



i"=i 



In rough terms, we can cut a tube in two at some geodesic time and find its amplitude 
by summing the products of the amplitudes of the pieces over all possible lengths of the 
intermediate loop (see Fig. 4). The motivation for the convention of taking the entrance 
loop to be unmarked and the exit loop to be marked lies in the fact that there are /" ways 
to glue the exit loop of one tube to the entrance loop of another — /" being the number of 
links on each loop — but we want to avoid factors of /" in the composition law. Suppose we 
mark the entrance loop of the rightmost tube in Fig. 4; such tubes would have amplitude 
^"^/",z'('^2)- We now can glue the two tubes together in such a way that the two marks are at 
any of /" positions relative to each other, so the amplitude for the resulting surface would be 
I" {Ni^iii[di)l" Niii i[d2)) ■ But that surface would still have the two marks on the intermediate 
loop, and we must delete them to get a surface of the type shown on the left side of Fig. 4. 
Thus we divide our last expression by /"^ to get Ni^i'/{di)Ni//^ii{d2) and sum over /" to get 



Since Eq. is just a matrix product, it is clear that Ni^i/(d) is the time evolution 
kernel of some Hamiltonian. We will find the continuum limit of this Hamiltonian and then 
calculate Ni^ir[d) by solving the corresponding Schrodinger equation. 

In terms of N{y,y',d) = Yji,u>iNi,i'{d)y\y'Y , Eq. (^) takes the form 



where the integral is taken along a contour around the origin. Making the change of variables 
X = |/~^ exp(— e^), we find 



Eq. (H). 



1 r dx 

N{y, y', di + d2) = —f —N{y, x, d,)N{l/x, y' , ^2) 



(33) 




- / diN{C.ON{-iX'.d) 

HI J-ioo 




N{C- a'e-fiC, r), d) = iV(C, C, d) - a'e^iC, r)j-N{C, d) , 
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where we have closed the contour to the left because the function N{—^, (, d) has a branch 
cut for positive real ^ > -y/r. Setting D = a'e'^d and taking e — 0, we find that A^(C, D) 
is a solution to the equation 

A^(C,/^) = -/(C,r)^V^(C,I^) (35) 

with the initial condition read off from Eq. (§T]), A^(C, C, 0) = 1/(C + C')-[1 In [0, this equation 
was derived for the case k = 2, but using only triangles in the discretizations, rather than 
squares; thus for this case we have a universality check. What is interesting is that the same 
equation holds for higher k, with the disk amplitude f{(, r) adjusted appropriately, and with 
a new scaling behavior for macroscopic geodesic distance D. The definition D = a'e^d is 
very interesting, and in a way unexpected, because it means that D has the dimension of L" 
instead of just L (recall a = /c — 3/2). 

In what follows it is convenient to introduce N{(,L',D), the inverse Laplace transform 
of N{(, D) with respect to C'- Since Eq. ( ^51) involves only C, and D, N{(, L', D) is 
also a solution of this differential equation, only with different initial conditions, namely 
A^(C,-Z^',0) = exp(— CL'). Our next step will be to find an approximate method for solving 
the differential equation and to use the solutions, as was done in |^, to determine universal 
functions describing the way random surfaces branch. These functions directly depend on 
N{L, L' , D), the inverse Laplace transform of N{(,L',D) with respect to (. N{L,L',D) is 
the continuum amplitude for tubes of length D with boundaries of length L and L', and it 
satisfies the initial condition N{L, L', 0) = 6{L — V). 



4 Calculating the tube amplitude 



A solution to Eq. ( pSf ) may be found by the method of characteristic curves: if the function 
Co(C,^) solves the ODE 

^Co(C, D) = -/(Co(C, D),t), Co(C, 0) = C (36) 



then ijj{(,D) = ^/'(^o(C) -D), 0) is a solution to Eq. (p5|). Geometrically speaking, ^ip^CD) is 
constant along a family of (non-intersecting) characteristic curves in C^-D space; Co(C) D) is 
the ^-coordinate of the point where the characteristic curve passing through ((^, D) meets 
the line D = Q. From Eq. (|36D it follows that Co(C) D) is implicitly determined by 

^^^77^ = ^. (37) 

Co f{i,r) 

Now A^(C, L' , D), the solution of Eq. (^) with initial conditions N{(, L', 0) = exp(— ^L'), is 
simply given by 



^In order to work with finite quantities in the continuum Umit, we rescale -/V(C, C'j ^) to absorb the factor 
1/e present in the lattice definition. 
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NiC,L',D)= expi-CoiC D)L') . (38) 

The inverse Laplace transform of this function gives the desired tube amphtude N[L, L' ,D). 
Thus, the problem is reduced to the explicit determination of (o{(,D). 

For k = 2, (q{(,D) can be expressed in terms of elementary functions. Indeed, after 
evaluating the integral in Eq. (|37| ) with f{(, r) = (2^ — y/r) \/ C + V^i and some algebra, we 
arrive at 



3 

Co = -Vr + -v^coth^ 



/3,M^ 1. VC + v^-Vl^' 



Unfortunately, it seems impossible to find the exact inverse Laplace transform of Eq. (|38|). 
Instead, we will look for the behavior of N{L, L', D) in the limit of a point-like entrance loop 
[L —>■ 0), as well as in the limit of a very large entrance loop [L oo). In the first case, it 
is sufficient to study the dominant behavior of Eq. (|39| ) for large and in the second case, 
for small C- 

Let us first consider the limit of a point-hke entrance loop. Substituting Eq. (|39|) into 
Eq. (pSf) and expanding for large C, we find 



+ o(i/C) 

(40) 



Performing the inverse Laplace transform, we find 

N{L 0,L',D) = _Le^'v^Ae-t^'^-^^(v^^^^^^) +0(L° 



(41) 



This formula is valid for L <^ D^, r~^/^. Nothing is assumed, however, about the relative 
magnitudes of and r"^/^. In the calculation was performed in the regime L -C <C 
r^^/^ so that the total area of the surface could be sent to infinity while D and L' were kept 
finite. Our Eq. (^if) generalizes the result of and reduces to it when expanded for small r. 
For A; > 3, however, we cannot find an explicit form of (o{(,D) analogous to Eq. (|39|) . For 
that reason we will work in the limit 1/C ^ D^^"' <^ r^^/^ (recall a = k — 3/2), and expand 
Co in powers of 1/C and r. 

Since /fc(^, t) = (—1)^ ^-^(l+0(r/^^)), as our first approximation we have from Eq. (|37|) , 

Co(C, D) = {{-1)'D + C^y'^^ + 0{r) . (42) 

From this formula we can see a sickness in the case of odd k: as D ^ Q^" from below, 
Co(C) D) — > oo. What is happening is that the solution to Eq. (^) with r = is running off 
to oo in finite time. Eq. (^) becomes complex for D > C"*^, but the real solution to Eq. ( p6D 
just stays at oo. Letting r become finite does not help the situation except when the initial 
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( is less than the rightmost zero of f{C,T), which is at C ~ V^- Since N{(,L',D) vanishes 
for large enough real (, the inverse Laplace transform N{L, L', D) cannot be positive definite 
for small L. This seems to lead to meaningless results, and we will not consider the limit 
of a point-like entrance loop in the odd k models. It may be that if we worked with fixed 
L from the beginning, instead of with Laplace-tranformed equations like Eq. (|35D, we could 
avoid the sickness we have observed. 

The even k models avoid this sickness, as is clear from Fig. 5 where we show the qualitative 
behavior of Co(C)-D) for A; = 3 and = 4 and finite r. To determine the baby universe 
distribution function in the case of even /c, we need to continue the expansion of Eq. (|42|) in 
powers of r and find the coefficient of the leading fractional power. First we expand 



1 



A; -3/2 „A; + 2n-3/2 , , ;,_,/23A; - 5/2 , ^, 



(43) 



where a„ and h are numerical coefficients which can be found from Eq. ([T8|) . Substituting 
Eq. ( ^31) into Eq. (^) and integrating term-by-term, we obtain the following relation: 



Co(C, D) ={D + Cn'"" + (D + C'^Y^-^I'^ (Co— 2 - C-'^-') + . . . + 0{r 



^-l-l/o■ /^-3cr-2 



-30--2 



0{ 



T 



(44) 



Now the expansion of Co in powers of r can be found iteratively. The coefficient of the 
leading fractional power, r'^'"^/^, is actually obtained after one iteration. Thus, the desired 
expansion has the form 



Co(C,^) = {D + CT'^'' + Tf 



[D + C 



{D + C 



-crN2+l/cr ^-3o--2 



- C'^-'' {D + C 



--y^-^''' + . . . + 0(r*^-i) 

+ 0{t^) . (45) 



Substituting this into Eq. ([38|), we have 



iV(C,L',D)=e 



(^ + C 



+ 



-0(r 



{D + C 



-(T\2+l/o- A-30--2 



- C-^--^ {D + c 



_<^N-l-l/(7 



+ 0(r* 



(46) 



Now we expand this for large C and inverse Laplace transform term-by-term, discarding the 
0{C^) pieces, which give zero length terms. We find 



N{L ^ 0,L',D) 



L 



L' 



Oir) + . . . + Oir^-^) + 



-fiDL'{{2k - 2)0^'" + L'y^-^l^ + 0{t'') ] + 0{L^''~^] 
where x = L' /D^^'^ and 71 = —hja^. 



(47) 
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5 The branching structure of random surfaces 



We now wish to extract from iV(L, L', D) a universal function which describes the branching 
structure of a large surface of planar topology and one boundary loop of length Lq. We will 
find that such surfaces have lots of little protuberances — hair, if you will. Let R{D) be the 
part of the surface which is at a geodesic distance D or less from the boundary loop. Let 
p{Lq,L,D) be the distribution of boundary loop lengths for R{D): on average, R{D) has 
p{Lq, L, D)dL boundary loops with length in the interval {L,L + dL). This quantity was 
introduced and calculated for = 2 and Lq — > in 0. The strategy for calculating it is as 
follows. The partition function for all surfaces is the disk amplitude F{Lq)/Lq. Let us think 
for a moment of the discrete version p(/o, /, d) of p{Lq, L, D). This function is the statistical 
average over all disks with boundary length Iq of the number of loops with / links at geodesic 
distance d from the boundary. According to [0 , in the limit of large disk area 

p(/o, /, d) = lim ^ y n . (48) 

(^j [GJlo] 

Nig^i{d)Gi/l generates only surfaces with at least one boundary loop of the desired sort. 
Moreover, if there are p such boundary loops on a particular surface, Ni^^i{d)Gi/l will generate 
that surface p times: each time, the tube generated by Ni^^i^d) will have a different one of 
its boundary loops left open for Gi/l to plug. Hence NiQ^i{d)Gi/l is the sum over surfaces 
generated by Gi^/Iq of the weight of each surface times the number of boundary loops at 
distance d and of length / on that surface. To understand the presence of the puncture 
operators d/dr, think for a moment of surfaces with fixed area instead of fixed r. To get to 
the fixed area representation, we would separately carry out inverse Laplace transforms on 
the numerator and denominator. The large area behavior of the fixed area quantities would 
be controlled by the leading singular term in the small r expansion of the fixed r quantities. 
Inserting enough puncture operators d/dr in the numerator and denominator of Eq. ( ^Sf ) 



to eliminate the leading analytic terms in r thus provides a convenient way to isolate the 
terms that survive in the large area limit. In order to obtain a generalization of Eq. ( ^HD 
to disks of finite area, we would have to replace the numerator and denominator by their 
inverse Laplace transforms. 

The continuum limit of Eq. (^) is obviously 

— , [N{Lo,L,D)F{L)/L] 
p(Lo, L, D) = Win ^ ^^/^xn . (49) 

Let us first consider the limit Lq — where the entrance loop is shrunk to a point. 
We expand 
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Fk{L,T)/L 



+ 0(t) + ... + O(r^-i) + 72L^-^/V'=-i/2 ^ 0/ 



(50) 



£^k+3/2 

where 72 is another numerical factor and we have adjusted Fk{L,T)/L by an overall multi- 
plicative factor for convenience. In this expansion, as in Eq. (0), only the 0(r°) term and 
the leading nonanalytic term are relevant to p{Lq, L,D). This is because both the numerator 
N{Lq, L, D)F{L) / L and the denominator F{Lq) / Lq of Eq. ( ^9|) have k leading terms analytic 
in r, followed by an 0{t''~^^'^) term. The leading terms are deleted by making k punctures. 
If one made fewer punctures, a positive overall power of Lq would make p{Lo,L,D) identi- 
cally 0. If one made more punctures, the leading nonanalytic term would still be dominant 
as r — > 0, so the final result would be unchanged. 

In the r — > 0, Lq — > limit, the only finite physical quantities in p(Lo, L, D) are L and D, 
so the only possible dimensionless scaling parameter is x = L/D^^". Thus p{Lq 0,L,D) 
is a function only of x, up to a dimensionful overall factor: 



p{Lo-^0,L,D) 



1 



7i 



72r(cr 



-X 



'<T-2 



(2(T + 1 + x) + 



r(a 



(51) 



where cr = A; — |, and the numbers 71 and 72 can be calculated explicitly: 



23-2fc 

" (6A:-5)(2A;-3) ' " (2A; - 1)!!(2A; - 3)!! ' ^^^^ 

The terms in square brackets in Eq. (^TJ) separate beautifully as A; ^ 00: the first term 
gives a non-integrable divergence at x = 0, and the second term gives a Poisson distribution, 
normalized to one and peaked at x = ex. 

The first term shows that the random surface has huge numbers of protuberances whose 
circumferential length is small compared to their geodesic distance from a given point; this 
is what we mean by the surface being hairy. The profusion of "microscopic" boundary loops 
suggests that as a one- dimensional universe propagates through geodesic time, it emits a 
divergent number of baby universes. 

For large x, the second term in square brackets is dominant, and it shows that there 
is exactly one "macroscopic" boundary loop to our region, and for large k its length is 
sharply peaked about aD^I" . Thinking again of a one- dimensional universe propagating 
through geodesic time, we would interpret the one macroscopic boundary loop as the parent 
universe, which survives the emission of its numerous baby universes. 

To determine at what x the first term becomes significant, we note that 



-log 
cr 



71 



72r(a 



-X 



-cr-2 



(2(T + 1 + x)e" 



;i + s + logs) 



(53) 



ascr = fc — 3/2^00 with s = x/a. This is positive for s < Sc and negative for s > Sc, where 
Sr satisfies 



15 



1 + Sc + log Sc = , 



(54) 



the solution to which is Sc ~ 0.2784645428. Let us define p(s) by 



p{s) ds = p{Lo 



0,L,D)dL 



■X 



cr-2 



{2a + l + x) + 



e dx . (55) 



72r((T) 



r(a + l) 



It follows from the above discussion that, as A; ^ oo, p{s) converges in the weak sense to a 
distribution which is +oo for s < Sc and 5{s — 1) for s > Sc- In Fig. 6 we show plots of p{s) 
for k = 2, where there is no separation of the "microscopic" and "macroscopic" terms; for 
k = 6, where the separation is significant; and for k = 100, where the convergence to the 
limiting case is very clear. 

Large k corresponds to the central charge going to — oo, which is held to be the semi- 
classical limit where quantum fiuctuations vanish and the surface is smooth, with constant 
scalar curvature. In the thermodynamic (infinite area) limit, the surface would be a plane, 
and we would have p{Lo —>■ 0, L, R) = 6{L — 2'kR), so p(s) = 5{s — 1) with s = L/ {2ttR). It 
is intriguing that this classical term is present in semiclassical limit we found for p{s). The 
different scaling law. 



and the profusion of microscopic boundary loops for s < Sc, are striking features of the 
quantum case which we cannot conceive of predicting by quasi-classical arguments. The 
scaling law Eq. (|56D might seem to be an artifact of the combinatorics: on the discrete 
surface, the matter fields are incorporated into the manifold by polygons with different 
numbers of sides, so the notion of defining geodesic distance as the minimal number of 
polygons one must traverse to get from point to point is suspect. So, couldn't we just define 
R = oD^I" j2Tx and claim that R is the "real" geodesic distance? The problem with this 
approach is that D enjoys a linearity property that seems essential to the notion of geodesic 
distance. Namely, if 7 is a loop each of whose points is a geodesic distance from a given 
point P on a random surface, and if 7' is a loop each of whose points is a geodesic distance 
D2 from 7, then each point on 7' will be a geodesic distance D\ +D2 from P. Any increasing 
function of D with the same property would have to be linear. The authors feel that a 
resolution of this question will have to come from a continuum formalism where the matter 
fields are more easily distinguished from the metric on the manifold. 

Another interesting limit where exact calculations are possible is that of an extremely 
long entrance loop, Lq —>■ 00. Here it is appropriate to calculate A^(C, L', D) in the limit where 
and r~^/^ are much larger than and L'. Now the characteristic curve equation, 
Eq. (^71), has a very simple approximate solution. 




(57) 
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so that, from Eq. ([38|), 



iV(C, L\ D) = l-CL' + L' Df kiC, t) + ... (58) 
Performing the inverse Laplace transform, we obtain 

N{L oo, L', D) = L'DFk{L, r) . (59) 
Using this and Eq. (|S^) in Eq. (^5]), we arrive at a remarkably simple formula: 

p{Lo oo, L, D) = LoDL-''-^/^ . (60) 

This distribution is valid whenever r^^/^ ^ Lq and both are much larger than L and D^^'^. 
The fact that the number of emitted baby universes scales as Lq could have been anticipated: 
this is simply due to the fact that a baby universe can split off anywhere along the parent. 
It is also clear why the number grows linearly with the elapsed geodesic time D, since in 
the discrete case the number of baby universes that split off at each step of evolution should 
be constant as long as the length of the parent universe does not change appreciably. Note 
that the distribution in L is again non-integrable for small L: the emitted baby universes 
are overwhelmingly likely to be microscopic. 

Another consequence of the preceding discussion is that 

N{L,L' ^0,D ^0) = L'DFk{L,T) . (61) 

This formula is interesting because it establishes a connection between the tube amplitude 
and the disk amplitude. N{L, L' ^ 0, D —>■ 0) is the sum over disks of boundary length L 
with a marked point located at a vanishing geodesic distance D ^ from the boundary. 
Thus, we expect that in the L' ^ 0, D ^ limit the tube amplitude reduces to the 
disk amplitude with a marked boundary point. This is indeed what happens, according to 
Eq. (|6TD . We speculate that Eq. ( |6ll) or some more general form of it could be used as the 
basis of a continuum derivation of Eq. (^5]). 



6 Discussion 



In ref. |TT[ a string field theory formalism for c = gravity was introduced. From this 
formalism Eq. (|35|) was elegantly derived. In fact, Eq. ( p5D was derived first in ref. via 
a careful combinatorial analysis of discretizations, and the string field theory was tailored 
to reproduce this result. In this paper we extended the combinatorial analysis to arbitrary 
discretizations with 2m-gons and established the validity of Eq. (^) for 2-dimensional gravity 
coupled to the (2, 2fc — 1) minimal models. This strongly suggests that the string field theory 
formalism of ref. encompasses all these theories. We choose a particular theory only 
through its disk amplitude /(C,'?"), which is the background value of the string field. 
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In solving Eq. ([35| ) we found a drastic difference between the even and odd k. In other 
calculations no such major differences were noted. For instance, the disk amplitudes of 
Eq. (^) are positive for all k. Thus, we may have the first indication of a serious difference 
between even and odd k occurring for spherical surfaces. It would be nice to understand a 
deeper reason behind this. 

The non-integrable divergence of p{L, D) for small L indicates that random surfaces are 
very hairy, even when restricted to spherical topology. In terms of the propagation through 
geodesic distance of a loop along its world-sheet, we take this to mean that tiny loops — 
the baby universes — are constantly splitting off the main loop. The conclusion that the 
microscopic baby universes are overwhelmingly more likely to split off than the macroscopic 
ones is quite intriguing in light of the large wormhole problem [|]. 

One might think of the microscopic baby universes as analogous to the soft photons which 
create the well-known infrared problem in the bremsstrahlung cross-section. Therefore, only 
the inclusive probabilities, where we sum over all possible splittings, are non-vanishing. We 
have verified that the probability for a loop to propagate any finite D along its world sheet 
without splitting is zero. 

The limit k oo (corresponding to the central charge decreasing without bound) bears 
a subtle relationship to classical gravity which we do not fully understand. Suggestions of a 
manifold that is smooth at length scales large compared to D^^"^ emerge from the k ^ oo 
limit of p{Lq — i> 0,L,D). But baby universes still play an important role in this limit, as 
evidenced by the divergence of p{Lq ^ 0, L, D) as A; ^ oo for L < ScCrD^^" . Perhaps the 
proliferation of baby universes up to this critical scale is related to the presence of many 
operators of negative dimension. 

It would be very interesting to study models with unitary conformal fields coupled to 
gravity, to see whether the surfaces they produce are more regular or more wild. To calculate 
the transfer matrix with the same methods as described here, however, one would need the 
disk amplitudes of the model with arbitrary boundary conditions on the matter fields, and 
these are not available even for such simple models as the Ising model. 
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Figures 

1. Small example of a one-step evolution of a loop on a random surface tiled with squares. 

2. Basic shapes for A; = 2. 

3. The first few terms contributing to one basic shape, drawn on both the polygonal and 
dual lattices. 

4. The composition law. 

5. Characteristic curves for /c = 3 (top left) and k — 4 (top right) above graphs of fsiCT) 
(bottom left) and f4^{(,T) (bottom right). 

6. The scaling function p{s) for /c = 2, /c = 6, and k = 100. 
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